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Abstract 

We study the spectral correspondence between a particular class of Schrodinger 
equations and super symmetric quantum integrable model (QIM). The latter, 
a quantized version of the Ablowitz-Kaupp-Newell-Segur (AKNS) hierarchy of 
nonlinear equations, corresponds to the thermodynamic limit of the Perk-Schultz 
lattice model. By analyzing the symmetries of the ordinary differential equation 
(ODE) in the complex plane, it is possible to obtain important objects in the 
quantum integrable model in exact form, under an exact spectral correspondence. 
In this manuscript our main interest lies on the set of nonlocal conserved inte- 
grals of motion associated to the integrable system and we provide a systematic 
method to compute their values evaluated on the vacuum state of the quantum 
field theory. 
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1 Introduction 



The existence of common symmetries underlying both ordinary differential equa- 
tions and quantum integrable models was first observed for the simple anharmonic 
oscillator Schrodinger potential in [TJ, from which functional relations characteristic 
to quantum integrable models were constructed. Soon after, in [2], the addition of a 
centrifugal term to the differential problem was shown to correspond to a twist in the 
boundary conditions of the lattice model. Moreover, this allowed to relate a differential 
eigenvalue problem to more general values of the Fock vacuum parameter in the QIM. 
The extension to excited states of the Virasoro module was achieved in [3] by intro- 
ducing a singular part to the potential with poles located away from the origin. The 
overall symmetry of the differential problem in all cases, although hidden, is encoded 
by the A\ = su(2) Lie algebra. The quantum field theory associated to this symmetry 
has been thoroughly studied in e.g. |IH6]. Besides the twisted XXZ spin chain, also the 
6-vertex lattice model, the massless twisted sine-Gordon integrable field theory and 
can be shown to be spectrally equivalent to this Schrodinger equation. 

The crucial ingredients in both frameworks are the Symanzik symmetry of rotations 
on the complex plane in the case of the ODE and the existence of a Yang-Baxter struc- 
ture for the QIM. These two elements, although not directly related, lead to equivalent 
functional relations - the Bethe Ansatz Equation, the T-Q or the Quantum Wronskian 
relations - which allow one to establish a parallel between the transfer matrix T(E) 
and Baxter's Q(E) operator to spectral determinants D(E) and Wronskians W(E) of 
solutions to the O.D.E., e.g., [?]. Consequently, the eigenvalues of Q(E) have zeros at 
the same points of the spectral parameter E as D(E). Thus, lattice models, interesting 
from the statistical viewpoint, and Integrable Quantum Field Theories (IQFT), possi- 
bly minimal models in Conformal Field Theories (CFT), can be conveniently mapped 
onto ordinary differential equations defined on the complex plane. 

The problem in discussion is not only closely related to physical applications but 
also has connections to more purely mathematical questions such as number theory. 
In fact, the eigenvalues of this problem can be combined in spectral Zeta functions, 
e.g. [8], allowing one to obtain some of the nonlocal charges [9]. With this in mind, the 
ODEs and integrable models can therefore be used to infer information about the Zeta 
functions. Another important byproduct of the investigation on the aforementioned 
correspondence came in [10J, with what is generally considered a reality proof for the 
spectra of a class of non-Hermitian PT-symmetric Hamiltonians [LTJ[T2]. Therefore, 
the study of spectral equivalences between ODEs and integrable models found fertile 
ground on the PT-symmetry community. This area received a great deal of atten- 
tion from a theoretical point of view and more applications of the so called IM/ODE 
correspondence to PT-symmetric Quantum Mechanics are natural. 

In order to investigate the connection of the integrals of motion in the QIM to 
the differential problem, a simpler extension consists of utilizing the supersymmetric 
= SM (2|1) Lie algebra. The advantage of this relies on the existence of a well 
understood quantum field theory, namely the AKNS soliton hierarchy [TSJCH]. In this 
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situation we have close links to the Integrable Perturbed Hairpin model and the Perk- 
Schultz model. Another positive aspect is that the sw(2|l) algebra, being one of the 
simplest supersymmetric extension of the Lie algebra, provides a good environment 
to test ideas in supersymmetry and superconformal field theory, which could play 
an important role in the fundamental descriptions of nature, according to AdS/CFT 
correspondence [T5]. 

In what follows, in order to fix the notation, we will present the class of ODEs of 
interest and review how to construct, from its symmetries, spectral functional rela- 
tions. With a few identification we can show that this expansion corresponds to that 
associated to the quantum AKNS field theory and that the coefficients in this series 
are the eigenvalues of its local integrals of motion I n . The nonlocal integrals of motion 
H n , usually more difficultly obtained, are constructed through a novel perturbative 
analysis of the ODE after a suitable change of variables and the first few are explicitly 
presented. 



2 The Schrodinger problem and functional relations 

Interested in studying an ordinary differential equation which could encode all the 
information about the quantum AKNS field theory, we turn our attentions for the 
moment, to solving the following differential equation, 

(~ + 4 (P 2 + P^y) + k, 2 (e'y + e 2 ^)) $(y) = 0, (2.1) 

defined originally on the whole real line — oo < x < oo for a restricted interval of 
the coupling < (3 2 < \. This is the equation presented in [13] if we identify the 
parameters there as p = 2tP,q = 2P,n = —2f3 2 . For our convenience we introduce a 
preliminary change of variables 

V = log {{-E)-^x M+l ) , <&(„) = (¥±^(3,), (2.2) 

and a consort reparameterization given by 

- 21+1 P=^, ^=^T. «=^- (2-3) 



4(M + 1)' 4(M + 1)' ' M + l' M+l 

This allows one to re-express the so-called Integrable Perturbed Hairpin model ODE 
2.ip in a more familiar fashion, namely, a generalized harmonic oscillator, 

d 2 + £_{£ + }) + x2M + axM _, _ E \ = ^ (2 4) 



dx 2 x 

defined on the half-line only, < x < oo, within the range M > 1, for which the 
harmonic oscillator M = 1 is a limiting case. The appearance of the a-term was shown 
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in [16] to correspond to breaking the su(2) symmetry into su(2\l), for reasons which 
will become clearer later. The complete knowledge about the spectrum {E^} can give 
information about the ground state of the associated quantum integral model. It is 
also possible to extend (\2A\i once more so that the excited states can also be accessed 
by the introduction of singularities at £ 2 ( M+1 ) = ^ in the potential [3l[T7] . 

The equation above is an extension of the x 2M potential, or of the PT-symmetric 
non-Hermitian (ix) 2M problem [TTJ[T2], so that when x — > ix the associated eigenfunc- 
tions and eigenvalues are given by ip{x) = ip(ix),\ = —E. The investigation of this 
family of potentials has generated a number of works on non-Hermitian Quantum Me- 
chanics and for the non-supersymmetric problem, a = 0, the reality proof for its spectra 
was given in [10] , whereas metrics and Hermitian counterparts for the — x 4 -potential is 
found in [IS]. More general non-Hermitian Hamiltonians have been studying outside 
the scope of the IM/ODE correspondence in e.g. [HH2T]. It is also interesting to note 
that when taking x — > ix, the part of the potential which has singularities away from 
the origin changes by moving the position of the poles, — > e ~ 2m ( M + l ) Zk , T n fa e 
event of M being a half-integer the position of the poles change sign, whereas for an 
integer value they remain unchanged, being still closely related to those of the original 
problem. 

The limiting M = 1 case usually requires special attention as some of the results no 
longer hold and have to be modified. However, the introduction of an a ^ term has 
a unifying character, as the energy can be interpreted to be rescaled by an ammount 
a, in this case. Therefore, the usual re-adaptation one has to do when evaluating the 
extended harmonic oscillator can be easily obtained by simply taking a = —E, 

ip{x;M = l,£,a -E,E = 0) = ip(x;M = 1, £,a,E), (2.5) 
ip(x; M = l,£,a = —E, E = 0) = ip(x; M = l,£,a = 0, E). 

This model (M = 1, a = 0) is sometimes denoted the Hooke atom, because it describes 
a particle in the presence of a centrifugal form alongside a linear Hooke force, is a bench- 
mark for being exactly solvable. In the integrable models framework, this corresponds 
to a coupling of (3 2 = |, known as the free fermion point because in this situation the 
corresponding sine-Gordon field theory can be mapped onto a non-interacting mas- 
sive Thirring model. The solution of the Schrodinger problem is given in terms of 
Laguerre polynomials and Hypergeometric functions, depending on the boundary con- 
ditions and eigenvalues, weighted by a decaying exponential factor. For general M/l 
the quantization problem cannot be solved exactly in closed form unless E = 0. In this 
case, 

ip(M,£,a,E = 0,x) oc 



i[)(M,£,a,E = 0,x) oc 



where L and U denote Laguerre and Hypergeometric functions. 



x i(i±(2m)) e -^rr L M+ 1+i 



M+l 21+1 / 2x + 



M+l±(2l+l)+q 
2(M+1) 



M + l 



M+l AI+1±(2£+1) /n M+l\ 

i(i±(2*+i)) e -W[/ — ] (2 6) 
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Still in the harmonic oscillator regime, we can make use of the fact that the differ- 
ential equation is also invariant under x — > ix, a — > —a, E — > —E, which can be 
seen as a PT-symmetrization of the problem. In fact, more generally, equation (12 .4p is 
invariant under discrete rotations by multiples of a certain angle on the complex plane, 
alongside a change in the constants a and E. The so-called Symanzik transformations 
for this system read [22] 

x -> u~ k x, a -»■ (-l) h a, E ->■ u 2k E, u = exp ( — — 

Since symmetries in general help one solving differential equations, we can use the this 
invariance to generate new solutions, suitably normalised, 

Mx;M,£,E,a) = (-l)* +1 l w^ 1 - e * ra V(w~ fc a;;M,^a; 2fc S,e ,fer a), (2.8) 




which solve the ODE for any fceZ. It is the unique subdominant solution, specified 
by the condition that it decays exponentially at infinity, in the region Sk, 



Sk 



arg(x) 



M + l 



< 



7T 



2(M 



(2.9) 



Once one knows a solution of the second order differential equation, the symmetry 
on the complex plane can be used to generate a new linearly independent solution, and 
one can use these two to write any solution of the problem. However, allowing (12 .4p 
in the complex plane requires a consistent choice of boundary conditions, as square- 
integrable wavefuntions can only exist in appropriate regions. The Stokes sectors in 
the complex plane define wedges where eigenf unctions exponentially decay and blow 
up in an alternating pattern. The problem we are mainly interested in throughout 
this manuscript is the radial problem, defined on the positive real axis and imposed to 
vanish at the origin and infinity, within the sector So- 

The k-th rotated solution, if)/- is subdominant in Sk and dominant in Sk+i- More- 
over, {ipk, V'fc+i} forms a basis of solution and W[ipk-i, V'fe+i] wm vanish if and only if 
there exists a solution decaying in both sectors. Thus, the solution ip = if) one may 
find for the problem is, for instance, linearly independent to ip±i, solutions on the com- 
plex plane rotated by an angle u ±l with respect to t/V In our problem, defined on the 
positive real line, with vanishing boundary conditions at the origin and at infinity, ipo is 
defined on the real axis, whereas ip±\ satisfy boundary conditions on certain regions of 
the complex plane. To be more precise, the appropriate boundary condition at infinity 
for our radial problem implies it decays exponentially according to 

/ X M+1 \ 

^(x; M, £, E, a) ~ M, £, E, a) = x~— exp I - — - J . (2.10) 

Note that when M is not an integer IN, the potential V(x) is not single valued, so a 
branch cut, chosen along the negative real axis, is necessary. 
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Another symmetry of the general ordinary differential equation fl2.4jl is that it is 
invariant under i — > — (£ + 1). Therefore, there are two linearly independent solutions, 

^ + (x;M,£,a) = ip{x;M,£,E,a), (2.11) 
if>-(x;M,E,a) = i/)(x; M, -(£ + 1), E, a). 

The parameter £ dictates the behaviour near the origin and the expressions above 
indicate that, although (I2.6P has a singular behaviour near the origin of the type 
x~ e , it can be easily analytically continued to a regular solution. Thus, by simply 
using this transformation one recovers square integrability, an important property when 
considering the Hermiticity of the differential problem. Near the origin the solutions 
behave like 

tp ± {x;M,£,E,a) ~ ^(x; M,£, E,a) = x^ ±2£ ^ , as x 0, (2.12) 

and one can use them as a basis of expansion for a solution in this region. Because 
the problem is not solved exactly in general for M ^ 1 one must instead resort to 
constructing solutions near the origin by the Frobenius method and using the WKB 
method to find asymptotic solutions which can furnish us with non-perturbative infor- 
mation about its analytical structure. 

The basis of solutions (p^ around the origin (I2.12p form a linearly independent set 
satisfying the following Wronskian identities, 

W[- +] = W[(p~(x; M, £, E, a), ip + (x; M, £, E, a)] = 2£ + l, (2.13) 

(2£ + l)ip(x;M,£,E,a) = W [0j+] {M, £, E, a)ip~(x; M, £, E, a) + 

+ ^[o,_](M,£,E,a)^ + (x;M,£,E,«), (2.14) 

with x-independent coefficients given by the Wronskians W[o,±] = W / [?/',V 9± ] of a so- 
lution tpo(x; M, £, E, a) with this pair, so that one can use them to expand a general 
solution in this region. The exact relation between them is given in terms of Wronskians 
and because ^[0,1] = W[-i,o] we can use a normalization condition of tpk so that 

VK[_i,i] (M, £, E, a)ip Q (x; M, £, E, a) = i[>-i(x; M, £, E, a) + ipi(x; M, £, E, a), (2.15) 

and in order to eliminate completely the x-dependence of this relation, one takes the 
Wronskian of this equation with the pair f ± , giving 

W hl)1] (M,£,E,a)W [0 , ±] (M,£,E,a) = u~ ! ^ ±2l W %±] {M , £, u~ 2 E, -a) + 

, a±(2i+l) , „ 

+ uj + — J — W [0 ,±](M,£,u; +2 E,-a). (2.16) 

One of the quantities in the left hand side involves the asymptotic behaviour in 
the sectors S±i of the complex plane and is therefore associated to the lateral problem 
and the "PT-symmetric differential equation; namely, VFr_x,i] vanishes if and only if, 
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for a certain value E n , the wave functions are such that ip±\ vanish asymptotically. 
Thus, this Wronskian encloses information regarding the quantization condition of 
this problem. On the other hand, the remaining object WV±] depends on the the 
asymptotic behaviour on the real line as well as near the origin, being associated to 
the radial problem we are interested in. 

Introducing 

C(E,a) = W hl , 1} (M,e,E,a), (2.17) 
D<*\E,a) = W [0i±] (M,£,E,a), 

the objects C(E, a), D(E, a) can be written as Wronskians of linearly independent so- 
lutions which are entire functions. Thus, these newly defined functions are also entire 
and we can use the Weierstrass factorization theorem to express the entire spectral de- 
terminants as an infinite product. A zero of D(E,a) indicates that the eigenproblem 
has a solution for which the eigenfunction corresponding to the eigenvalue E satisfies 
the correct boundary conditions. Therefore, we can conclude that this coincides with 
the spectral determinant, up to a normalizing factor, vanishing at the possible eigen- 
values of the differential equation. What is more, the definitions above furnish us with 
the analogue of Baxter T-Q relations (also in j9] and reference [1] therein) 

C(E,a)D^(E,a) = (u +2 E, -a) + tJ-^D^ {u~ 2 E, -a) , (2.18) 

whose zeroes E = Ej are characterised by 



(u- 2 Ef\-a 



£>(±) I co+ 2 E) 



(a) 



-a 



-J 2e+1+a \ (2.19) 



Part of the solutions for this equation correspond to C(E,a) and the remaining 
ones are associated to zeros of D(E, a). Because the problem defined on the real line 
is Hermitian the zeros of D(E,a) should be real. Therefore, any complex eigenvalue 
should be a zero of C(E, a) but there is no guarantee that the zeros of C(E, a) cannot 
be real. 

In order to use the Weierstrass factorization, we must take into account that a 
simple WKB approximation tells us that the zeros accumulate towards infinity along 

2M 

the positive real axis according to E n ~ n M + 1 , so if M > 1 we can establish a convergent 
product Z)( _ ) = D in the form 

°° ( E \ 

D(E, a )=D(0,a)H[l--^), (2.20) 
k =o \ E k J 

which leads immediately to an expression of the form of the Bethe Ansatz Equations, 

'Ei a) -u- 2 E { - a) 



fc=0 v Er - u*E) 
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In this case, they correspond to two sets of coupled equations, both for the values 
of +a and for —a. The problems defined for +a and —a have different properties, 
reflecting on the eigenvalues, or Bethe roots. In the differential equation framework, 
for one case we have a convex potential whereas for the other a double well potential, 
respectively [23J. For M — 1, we cannot guarantee the convergence of the infinite 
product above and changes are expected to be necessary. The relations (12. 5p . though, 
make it simple to obtain the correct equations for the T-Q and Quantum Wronskian 
relations. One should note that, rather than the Q(E, a)-operator itself, the spectral 
determinant D(E,a) should be identified with the quantity A(E,a), defined by 



Q (u ±2 E, -a) = u ±2i± ^A (u ±2 E, -a) . (2.22) 



(M±I)Tm+T _^M^ _ respectively; gives 



To obtain the correct normalisation factor to the Hadamard factorization we ob- 
serve that the solution of the differential equation at vanishing value of energy can 
be found exactly (12. 6p and by normalising (I2.10p and (12.121) for instance with ^= and 

— i — 21+1 r (-r 21 + 1 s 

--(21- 
2(M+1) 

D(0,c*)~(2l + 1)(^±±) ? N . (2.23) 

1 ' 1 ' M+2l+a+2 N 

2(M+1) 

and the expansion of the spectral determinant above (I2.20p has the expected leading 
behaviour 

2M 

log £>(£,<*) TVT -T7 , 2V^4^K , (2-24) 



cos 



(—) \ v ' \ r 



2M+1 ' 
2M 



well known in the literature. An important result of the spectral determinant expansion 
above is the appearance of the constants which are coefficient of the higher order 
terms in this expansion in powers of the spectral parameter. These correspond to the 
eigenvalues of the commuting set of local conserved operators in the quantum AKNS 
model. The reformulation here presented is well known now and allows one to transform 
our original differential problems into a set of coupled nonlinear equations. In the next 
section we will show how the coefficients in the expansion of the spectral determinant 
can be identified as the integrals of motion corresponding to the quantum integrable 
model. 



3 Spectral Determinants and Integrals of Motion 

It is an important result that one can transform our initial eigenvalue problem into 
algebraic equations but the simplification is only relative as the number of equations 
to solve becomes infinite. In any case, this alternative method which is based not 
on the differential equation directly but rather in the functional relations one can 
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construct from it allows one to extend the toolbox to solve the original eigenvalue 
problem as well as giving access to important information in an exact analytic form. In 
the thermodynamic limit, this approach in essence transforms the differential problem 
in nonlinear integral equations, usually referred to as NLIEs. This has a central role in 
our strategy to compute nonlocal conserved quantities for the QIM from an analysis 
of the ODE. 

The starting point is to consider the C-D relation f |2. 18j) . which is equivalent to the 
T-Q relation 

T(E, a)Q(E, a) = Q (co 2 E, -a) + Q (tu~ 2 E, -a) , (3.25) 

if we identify the spectral determinant D(E,a) with the quantity A(E,a) in ( 12.22ft . 
One then defines a new object, 

a(E,a)-u r, (3.2b) 

D(uj *E[a),— a) 

so that for the zeros E k (a) of A(E(a), —a) the associated Bethe equation may be recast 
in the form of 

a(E k , a) + 1 = 0. (3.27) 

This is actually the same system of two coupled equations, one for +a and one 
for —a. We make the assumption that T(E,a) and A(E,a), or equivalently C(E,a) 
and D(E, a) respectvely, are entire functions so that we can factorize them in terms of 
their zeros in a convergent product through a Weierstrass decomposition. Form such 
a factorisation we can observe that a (E, —a) a(E, —a)* = 1, if the energy levels are 
real. 

One can write the quantisation condition as 

oo / i ,J2 E(a) \ 

log a(E, a) = log (co 2 ^) + £ log g|g = 2m (m k + - ), (3.28) 

k=0 V 1 U E k (-a)J V ' 

so that for each eigenvalues E k there are infinite possibilities for the integer parameter 
mfc. When we take the logarithm of the expression above we introduce an arbitrariness 
in the choice of the branch, or the Riemann sheet, denoted by an integer m&, in such a 
way that eigenvalue E k is associated to one of the integers m k . A certain eigenstate of 
the Q operator is characterized by a unique set of integers {m k }. However, not every 
set of integers will correspond to an eigenstate of Q. 

We can replace the infinite sum over zeroes by a contour integral. The total contour 
can be split in one part containing the positive real line and one contour for the 
remaining complex plane. Using the reparameterization 



2M 



E = e^+T, (3.29) 
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the logarithm expansion of the quantity a(9, a) can be expanded as an integral which 
can be written as a sum of infinitely many residues plus a diverging contribution at the 
large circle [23JI2I3, after the inclusion of a term responsible for the leading behaviour 
following ( gjMJ) , 

loga(M) = log (co 2e+1+a ) - 4£re e + (3.30) 

V 2 M 

+ / K + (0 -w)log(a(w,a) + l)dw - / K + (6 - w) log ( — - — r + l) dw + 



dw, 



Cl Jc 2 



a(w, a) 

+ / K(6-w)log(a(w,-a) + l)dw- [ KJO - w) log ( — — + 1 ) 

J a Jc 2 \a{w } -a) J 



where the contours c\, c<i are placed just above and below the real axis and the exact 
forms of the kernels K±{9) were given in integral form in [SSIEZ]- These integrals can 
be written as infinite sums, so the kernel may be expressed as 

K±{0) = \{<p{e)± X {&)), (3.31) 

with 

tptf) = y j£ e -MMn tan(7rMn) _ y l e g(l-2n) CQt , (3.32) 

n=l n=l v ' 

xw = E^ 9M(1 - 2 " ) «' t (f M ( 2 "- 1 ))-E^ 2, " t ' m (i). < 3 ' 33 ) 

n=l n=l 

Therefore, the sums over the residues only give an asymptotic series expansion 
and from them we can construct the conserved charges. The NLIE for the KdV su(2) 
problem, with vanishing a, involves a function which coincides with the solition-solition 
scattering amplitude for the sine-Gordon model. The same NLIE was earlier derived 
for the Bethe ansatz equation associated with the XXZ model and by taking the ap- 
propriate continuous limit leads to the sine-Gordon field theory. Thus, the integrals 
of motion of the aforementioned models are related to each other. In this case, the 
scatering matrices are the ones that should appear in the quantum AKNS field theory 
and the conserved charges are equivalent in both models. 

After some identifications we can relate the integrals of motion to the integrals 
appearing in the expansion. In the case of the nonlocal ones, we have 



1 (M+l\ 2n M 



hn --^\~^) cos(TrMn) X (3 ' 34) 



( f e 2nMiw log (a(a, w) + l)dw- [ e 2nMiw log ( , 1 r + l) dw 
\J C1 Jc 2 \a{a,w) J 

' 1 

a(— a, w) 



+ / e 2nMiw log (a(-a, w) + 1) dw - I e 2nMiw log ( — + 1 ) dw ) , 



' Cl J c 2 
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hn-i = -- Mcos(— ^ ^jx (3.35) 

f e (2n-l)Mi W log ,/ w \ + l\ dw _ f e (2n-l)Mi W j f 1 + -A rfw 

y ci i C2 v. a (a,^) / 

- / e ( 2 «-D A «- log ( a (-a, w) + 1) + / e i2n ~ 1)mw log f — r + l) dw] . 

Jci Jc 2 \a(-a,w) J J 

and similar expressions for the local ones [TTJISH]. We notice that 

In order to have a consistent notation with respect to the literature it is convenient 
to make the following identifications H n = I<m and H 2n-i = I% n -\ for when a^O. For 
vanishing a, the contributions of the integrals coincide so the last two terms in the 
series above cease to exist and one recovers the results of [H46|[28]. A WKB analysis 
of (12.41) leads to an expansion of the spectral determinant in terms of local conserved 
charges which can be found in e. g |3][28] when a = 0. The local charges for the 
supersymmetric case are investigated in detail in [T7] but also have both odd and even 
contributions, the latter vanishing with a, as expected. 

The expressions above though provide a recursive scheme which is efficient to make 
numerical comparisons. It is also an important tool from a practical viewpoint, due 
to the difficulty of compute the objects I n and -especially- H n , the local and nonlocal 
conserved charges of the model. Having in mind the equivalence between Baxter's 
Q(E) operator and the spectral determinant D(E) we will show in the next section 
how the nonlocal integrals of motion can be determined in a reasonable way. 



4 Nonlocal Conserved Charges from the ODE 

We have seen that the supersymmetric model here studied has an important feature, 
namely, it interpolates the properties of the quantum KdV field theory or, equivalently, 
that of the sine-Gordon model and the XXZ spin lattice, which are all well understood. 
Besides being characterized by local integrals of motion of the kind previously presented 
for the AKNS analogue, there is another set of important commuting operators, the 
denoted nonlocal integrals of motion. Originally defined to be the coefficients in a 
Taylor expansion of the transfer matrix operator, e.g. [5||29]. 

oo 

T(A) = 2 cos(2ttp) + G n \ 2n , (4.36) 

n=l 

with A being the spectral parameter related to E according to E — X 2 p, where p = 

2^(M + i)imryi T ) 2 . 

An alternative way of expressing the nonlocal conserved charges, in the case of an 
underlying su(2) symmetry, consists of defining them as H n in terms of the coefficients 
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r(i-/3 2 ) 

of an expansion of the Q-operator in powers of the variable A m — -■ The dual nonlocal 



integrals of motion, specified by H n ((3 2 ) = H n ( J , give rise to the following expansion 

of the operator A(E), which corresponds to the spectral determinant D(E) of the 
corresponding Schrodinger problem. This means we can write down an asymptotic 
expansion, possibly with limited radius of convergence, 

togD(^~£(jjrrr) U-Ey n{M+l \ (4.37) 

n=l ^ ' 

for large values of the spectral parameter E. 

However, we have seen that the addition of the a-term brought another set of local 
integrals of motion, characterized by even indices, to the pre-existing ones, specified 
by odd indices only. Therefore we would expect the same to happen with the nonlocal 
charges. This implies having an expansion 

00 

T{\,ot) = T(0,a) + ^G f A n , (4.38) 

n=l 

which reduces to ( I4.36P when all even powers of A vanish with a. So the first nonlocal 
conserved charge Gi appears in power of A, rather than G\ at order A 2 , and has no 

2 

equivalent when a = 0, unlike G\. Moreover, the dual nonlocal charges also come as 
twice as many, 

log D(E, a) ( WTl) ^(-E)- n ^, (4.39) 

n=l ^ ' 

and the first integral motion in the expansion would have no counterpart in the su(2) 
scenario and Hi would be the first nonvanishing contribution since no Hi would be 
present in that case. 

The nonlocal integral equations allow us to compute the charges numerically but 
analytic expressions are difficult to obtain, particularly in the case of the nonlocal 
integrals of motion. In [IH6] the G n are formally computed in generality for a = 
and exactly evaluated for first two contributions in the expansion. It would be in 
principle advantageous if one could determine these important objects by investigating 
the associated ODE with nonvanishing a by the existing spectral correspondence. In 
fact, here we show that these properties of the operators A(E, a) and Q(E, a) can be 
obtained in a systematic fashion by investigating the same original differential equation 
(12 .4p . Although the construction of the local integrals of motion can done directly from 
the latter, in the case of nonlocal charges some reparameterizations prove to be very 
convenient. For that, we introduce the new variables, 

log (±e~^x) , E(s) = -^=ip(x), (4.40) 

'x 



where the e is related to the spectral parameter through 

T = e M = (±E)- M P. (4.41) 
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This transforms the problem into an equivalent one, but free from singularities at the 
origin again, 

d 2 /2£ + l^ 2 



ds< 



E(s) = 0, (4.42) 



which corresponds to a potential expanded in terms of powers of J 7 . This can be 
associated to the AKNS problem by taking M = ^f-,^ = 2-p — |, a — 4-^r. 

For simplicity we can work with the choice e = E~^A~\ otherwise, we must replace 
s by s + z|, closely related to the associated "PT-symmetric problem. This struc- 
ture tells us that for small values of J 7 , i.e., large values of E, this equation can be 
solved perturbatively, order by order. The unperturbed potential, with T = 0, which 
also corresponds to s — oo, has two linearly independent solutions given by Bessel 
functions, 

/ of + 1 \ 

E±(s) = A4 T M ± J ± up (e s ) , (4.43) 

satisfying required boundary conditions of vanishing at s = ±oo and s = oo, respec- 
tively, as well as the familiar relations, with Af± = 7r ± 2r (^j-^ + ^") T > 

W[E_( S ),E + ( S )] = 2, 

E(s) = l -D + {E,a)H + {s) + ±D-(E,a)Z-(s). (4.44) 

In fact, the solutions E±(s) have first order contribution of the type v 9± ( a; ) i n fl2.12j) . 
after the appropriate change of variables. 

As s — > — oo, we have that E + (s) — > 0, so we can establish the proportionality of 
the spectral determinant D(E,a) = D(~'(E,a) = WV+i {E, a) with the wavefunction 
in this limit. However, expression f |4.44j) is generally true in this regime and does 
not rely on any interpretation of the D ± (E,a) as spectral determinants which vanish 
for appropriate boundary conditions. In fact, in the present context the transformed 
differential equation does not exactly constitute an eigenvalue problem, albeit closely 
related to one. Therefore, these Wronskians are simply the projection of a solution in 
the aforementioned basis and the problem of determining them relies on a solution of 
the differential equation. 

Because the latter can be seen as a perturbed equation in powers of the parameter 
J- we can assume a general solution of the form 

E(s) = E (+) ( S ) + E ( - ) (s), (4.45) 

oo 

E(±)( S ) = A<t\s) + £j*Af \s), (4.46) 

3=1 

where in lowest order Aq (s) = E±(s) is the solution of the unperturbed equation 
(J 7 = 0) satisfying the necessary boundary conditions. Besides, we have added all 
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the remaining perturbative contributions so that the result is expressible as a series 
in integer powers of J 7 . Restricting for now our analysis to the ground state of the 
Q-operator for the perturbed Hamiltonian ( I4.42p we can expand the wavefunction in 
powers of T and solve the associated Schodinger equation perturbatively at each order, 
first, 



(4.47) 



r h + sm\ 2 r s 
A?\s) = =F a S T 00 jf e^yj ±2 _^(eyy d y + 

then, by solving the resulting equation as the coefficient of J-" 2 , 

A?\s) = =F -4zTl Et(s) / e2(M+1)yj ± ^( e2, ) 2 ^+ ( 4 - 48 ) 
- f r (±^) M*) jT e^A 1 {y) J ±2 ^ {?) dy + 



1 7T 

± 



2 COs(7T^ 



£±(s) / e 2(M+1)?/ J_2^i (e*) J^±i (e y ) + 
E±(s) f e( M+ *A 1 ( l /)J T s i (^) 



2 V 2 

and so on, successively, for higher orders. We see that A 2 (s) depends explicitly on 
Ai(s) and the same pattern occurs at all orders: to determine the solution at each 
order it is necessary to know the solutions at previous ones. In fact the solution for a 
generic order can be found in closed formula described by the expression, 

r (M±l\ r s / x 

Af)( S ) = — 1 ^ 1 ^-(s)J s (e«0 (e^ASd,) + oAg(y)) dy 

The constants of integration are appropriately fixed to satisfy the required boundary 
conditions so that contributions which do not behave appropriately are excluded form 
A fc (s). 

When a vanishes, the contribution of A^s) is reabsorbed in the unperturbed solu- 
tion. Curiously, the contribution A 2 (s; a = 0) is almost the same as A 1 (s; a = 1) and 
that is because V ~ a e < - M+1 ) s J r + e 2 ( A/+1 ) s J-' 2 , i.e. when the supersymmetric param- 
eter vanishes the first contribution, coming at the order J 72 , is the square of the first 
contribution when such parameter is nonzero, at first order in T . 
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The expressions above seem cumbersome but from them we can see that it is con- 
venient to introduce the following two quantities, 



My) 



2 F 3 (£ + 1, f (M + 1) + g^; g±l + 1, f (M + 1) + + 1, 2(1 + 1); -e 2 ^) 

2 2£+l e -2((Af+l)n+2^+l) 2/ (( M + ^ + 2 £ + l)f (^±1 + l) 2 

KM _ ^ (j 1 (M + 1); f (M + 1) + 1, 1 - U±L, 1 + -e 2 ^) w 

nlyJ " (M + i)nr(i-^i)r(i + ^i) ' 1 } 



which arise as the indefinite integral forms involving Ja+i (e y ) 2 and Ju+i (e y ) J 21+1 (e y ) 

2 22 

in the expressions for A n (s) above. The equivalent result for J_«+i (e y ) can be triv- 
ially obtained by just replacing £ — > —(£-\-l) in the former. Taking the lower integration 
limit to be Xi = 00, or S{ = 00, and evaluating it near the origin, 1 = or s = - 00, as 
this is the region we are interested in the wavefunction, we obtain the quantities 

r = r(i-f(M + i))r(^i + f(M + i)) 

2 v ^r(i-f(M + i))r( 2 fi-f(M + i) + i)' 1 ' 

K = r(f(M + i))r(i-f(M + i)) 

2 v ^r(i- 2 fi-f(M + i))r(i + 2 fi-f(M + i))- { ' 

These are related to each other in a very simple way, 

J n cosf ((M + l)ra-l) 



K n cosf ((M + l)n-l + (2£+l))' 
and can be easily related to the quantity I2 = H\ calculated in [5], 



2 / 2 M ^ 2 



j2 = ; sinM,r i^TTj (452) 

with the sin contribution arriving because we are expanding in T = E i~ . These 
constants play an important role in the theory as they will then appear when one is 
evaluating the solutions at each order. 

We have seen that in every order there are contributions in the direction of both 
£±(s) and one can collect them and determine an expansion for the then called spectral 
determinant. This implies that, up to the Weierstrass normalisation factor which fixes 
D(0, a), an expansion of the spectral determinant in the form 

D(£,«)~l + yF n J n , where F n = lim " / ; . (4.53) 

£1 ' AT 00 



In this limit we are interested in, the above expressions allow us to compute the 
quantity Fi, appearing as the coefficient in the first order expansion in terms of E ~ 



14 



and it turns out to be given by a particularly simple expression which is proportional 
to the constant J\ we have just defined, 

Fx = a-J x . (4.54) 

This object has a simple form and does not contribute in the su(2) limit because it is 
proportional to the supersymmetric parameter a. 

Solving the second order term is trickier as it depends on the solution Ai(s), which 
is to be integrated in the object 

Bi(s) = [ e<- M+x)v &x(y)J ± up (e y ) dy. (4.55) 

The evaluation in this case, when lim^.oo, show that the factor F 2 , although having 
a more complicated form, once more the solution can be written down as a function of 
the fundamental quantities J n , or K n . It yields that 

F2 - ~J2 + a A7T 2 (2e + 1) (P l + b 2 ), (4.56) 

where the remaining constants 6j, which arise from the evaluation of Bi(s) at — oo, 
involve functions of the kind 5F4, 

= 2T(-M)T(M + i + l) 

my (~ — m) r(— m +£ + 1) ' v ; 
7rr(M + i)r(-M-f-i) F 

A e {M + 2£ + 2) 2 r (£ + |) 2 T(-M -11- 2)Y (-M - £ - \) 5 ^ 
and we have used for short the quantities 5F41 and 5-F4 2, respectively as 

5 F 4 (^,-M,-^-,-£,£+l;^-M,l-^-,-M-£,-M + £ + l;l 

5F4 [m + 1, 1 + 1, y + I + 1, M + i + ^, M + 21 + 2; £ + ~ ^ + £ + 2, M + £ + 2, 2£ + 2; 1 

The appearance of these hypergeometric functions already in ifi make evident the 
complications presented by the introduction of the supersymmetric term since for a = 
such contributions would appear at higher orders, starting at H 2 , or G 2 as in [9]. The 
present results, though, do not rely on any knowledge concerning the properties of Zeta 
functions. 

Now we are finally in a position to evaluate the first few terms in the expansion 
coming from (14.531) . so that the logarithm expansion of the spectral determinant takes 
the form 

log D(E, a) ~ FF X + T 2 (f 2 - ^pj + (4.58) 
(f 3 + 3- - F 2 F^ + T 4 (f 4 - S _ F 3 Fx + F 2 Fl 
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In the limit of vanishing a, the odd (or half integer, according to other notations) 
components vanish, as expected from the sw(2)-symmetric theory, for which only even 
(or integers, alternatively) terms are present. For this situation, the first contribution 
is given simply by 

/ nAI \ 2 

^ = (-1)- (M+1) (tttt] 4 (4.59) 



M + l 



with the standard identification 



/, = ft = (M + im -(2M + l ^'^ffi . (4.60) 

in agreement with f)4.39p . In possession of the expressions (I4.53P we can keep on 
with (I4.49p and determine higher contributions to the solution expansion and therefore 
extract higher nonlocal integrals of motion. Those terms however become cumbersome 
with simplifications when a = 0, as even F n disappear form (14.581) . 

Nevertheless, it would be interesting to verify if the method currently under investi- 
gation provides any information about the integrals of motion evaluated away from the 
vacuum of the Q operator. It was shown in [3] that higher level eigenvalues of the Bax- 
ter operator can be obtained by adding to the potential poles located at x 2( - M+1 ^ = Zk- 
For the first excited state, for instance, 

(2£ + l) 2 - 4M 2 
« = i '— . (4.61) 

The corresponding part that needs to be added to ( I4.42[) is also expressed as a power 
series in the parameter J 7 , 

U L (s) = 4(M + 1) X: E mM + , 1} " l) ^^J^^ , (4.62) 

k=l j=l Z k 

but the expansion is in rational, not integer, powers - unless M — 1. Therefore one 
needs to accommodate this feature in the solution series according to 

oo oo 

E( S ) = A (s) + J2 jrm ^ T ^> ( 4 - 63 ) 

3=1 3=1 

where Tj are the contributions due to the presence of the higher level potential (L ^ 
0), involving non- integer powers of the perturbation parameter. Here we simplify 
the notation by eliminating the upper indices ^ and present only the positive ones, 
reminding the negative ones are related to these by a simple change in the sign of the 
combination 2 ^-. 

In principle, there are infinitely many Tj because this is the expansion of the Ul 
potential and investigating the behaviour of the this introduction, in a similar way as 
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we have just done for L = 0, consists of finding solutions to the Tj(s) in the expansion 
(I4.63p . It turns out that these terms can be solved in closed form for a generic j-th 
coefficient, 



TjOO = S_(s)(2j(M + l)-l)(M + l)(2£+l)r 



2£ + r 2 



X 



x 



^— ' ~J / 2 2 

fc=l ^fc ■ / *J2 



Once more, in order to determine the contribution of the excited states of level L 
to spectral determinant expansion we must evaluate the wavefunction in limit where 
we take s — > — oo. In this case, the integrals, which have a similar form as those in 
(I4.49p and (14.491) . can be evaluated. With this expression at hand we have determined 
how the introduction of the poles Zk in the ordinary differential equation affects the 
spectral determinant expansion, 

t,(.) = r (i - *±i W («•) ^ M ±^ (m i)r (*±1) ' x 



1 r (\ (3 - 2j(M + 1))) r (|(2j(M + 1) - 3) + £ + l) ^ 1 
20Fr (|(3 - 2j(M + 1)) + 1) r (|(3 - 2j(M + 1)) + £ + l) ^ 



V 2 / 2 4 cos(vr£) 

1 T (i (3 - 2j(M + 1))) r (|(2j(M + 1) - 2)) " 1 
2 v ^r(i(3-2j(M + l))-£)r(i(3-2j(M+l))+£+l)^4' 1 " J 

evaluated at s = — oo. Therefore, the expansion suggested allows us to know the exact 
form of each of the higher level coefficients in the expansion. The complete evaluation 
of the function above requires, though, knowledge about the exact positions of the 
poles, which is known exactly for L — 1. However, even if we include the expression 
for z\ the overall contribution of the Tj does not appear in D(E,a) in the precise 
form one is interested, namely integer powers of J 7 . This prevents us from determining 
the effect of excited states on the nonlocal integrals of motion, indicating the above 
parameterization is not good to investigate the effects of higher levels L > 0. 

The situation is different when M — 1, with the introduction of poles contribute to 
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lower orders according to 

U L (s) = - e 2s F + ^ e " Sjr3 + ■ ■ ■ > ( 4 - 66 ) 

fc=l Zk k=l Z k 

so that the lowest contribution comes simply as a shift in the term of order J 7 in (14.42)) . 
In this situation, it is enough to make the following replacement in F\ and F 2 , (14.541) 
and (j4.56p . respectively, 

L OA 

a^J2-, (4.67) 

which for the first excited state can be determined analytically with ( I4.61|) . If the 
location of the poles Zk(a) are known for when q / 0, then in order to evaluate the 
first contribution of the excited states one needs to take a — > a + Vf , -^-r instead. 

A^fc=l z k (a) 

The fact that the poles are located away from the origin guarantees a well behaved 
modification to the supersymmetric parameter. Higher contributions can also be calcu- 
lated but to the best of our knowledge these objects have not been explicitly presented 
for a direct comparison. In fact the study of excited states deserves separate attention. 
The appearance of the new terms to the original potential at orders J-" 3 , J 75 , ■ • • does 
not allow us to extract the higher level contributions from the vacuum in a simple way. 
The solution for these cases is obtained by identifying Tj(s) with Aj(s) and following 
the same procedure for the modified equation at each order. 



5 Conclusions and Perspectives 

The spectral equivalence between sw(2|l) supersymmetric spectral determinants 
D(\) associated to a class of ordinary differential equations and the eigenvalues of the 
Baxter Q(A)-operator of the quantum AKNS field theory has been investigated. We 
have presented a framework capable of proving us with the nonlocal conserved charges 
associated to the quantum model, H n in order to determine the whole NLIE expansion 
together with the local conserved quantities I n . This constitutes an alternative and 
efficient method which allows one to determine both odd and even charges, which are 
missed out from the spectral Zeta function approach of [9], and extend the results 
in [13, 28J for a supersymmetric potential in an elegant way which can be reduced to 
the previous case when the deformation parameter, responsible for taking an su(2)- 
symmetric problem into an su(2|l) analogue, is set to vanish. 

Many questions remain open, especially regarding the understanding of the higher 
level IM / ODE correspondence for algebras beyond the local integrals of motion associ- 
ated to su(2) and su(2|l) [3|[TT] . Here a first step is taken in order to also accommodate 
the nonlocal conserved quantities. The introduction of poles in the complex plane may 
destroy in some cases the reality of the spectrum for a certain ordinary differential 
equation. Therefore, the excitation of higher Q(A)-states may be seen as a mechanism 
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of breaking of PT-symmetry for a non-Hermitian Hamiltonian, a phenomenon which 
deserves attention on its own. Once we have established this equivalence between cer- 
tain two systems also for the excited states, information extracted in one side may be 
valuable on the other. 
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